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Abstract 

(N 

The role of the reaction e + e~ — > e + e~ir + ir~ in the pion form factor measurements 
via radiative return method without photon tagging is studied in detail. The analysis 
is based on the developed Monte Carlo program EKHARA, which ingredients are 
also presented. 
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1 Introduction 



Radiative return method of the hadronic cross section extraction from the 
measurement of the cross section of the reaction e + e~ — * hadrons + photon(s), 
proposed already some time ago [1], is currently being used by KLOE [2] and 
BaBar [3] providing very precise experimental data. Further improvement in 
accuracy is crucial for predictions of the hadronic contributions to a M , the 
anomalous magnetic moment of the muon, as the error on the hadronic con- 
tributions to a M may obscure possible new physics signal, seen as a deviation 
from the Standard Model (SM) predictions. The same information is essen- 
tial for the evaluation of the running of the electromagnetic coupling {aQED) 
from its value at low energy up to Mz as the present error on the hadronic 
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contributions is too big to fully profit from the data of the future ILC (In- 
ternational Linear Collider) running in the gigaZ mode. For recent reviews of 
these subjects look [4,5,6]. 

The extraction of the cross section a(e + e~ — > hadrons) from the measured 
cross section cr(e + e~ — > hadrons + photons) relies on the factorization 

da(e + e~ — > hadrons + 717) = H da(e + e~ — > hadrons) , (1) 

valid at any order for photons emitted from initial leptons, where the function 
H contains QED radiative corrections. This function is known analytically, if 
no cuts are imposed, at next to leading order (NLO) and has to be provided 
in form of an event generator [7,8] of the reaction e + e~ — > hadrons + photons 
for a realistic experimental setup. 

Let us focus on the most important process, where the 'hadrons' means just 
7r + 7r~ pair. This reaction gives the dominant contribution to the hadronic 
part of a M as well as to its error. In case the photon(s) are not measured and 
only charged pions are tagged, there exists a number of possible backgrounds. 
It was pointed out in [9], basing on integrated over the whole phase space 
analytical formulae and containing contributions from diagrams (a) and (d) 
in Fig.l, that the reaction e + e~ — > e + e~7r + 7r~ can give sizable contributions 
to the radiative return process, especially for low invariant masses of two pion 
system. To examine this contribution for a realistic experimental setup, a 
Monte Carlo program EKHARA was developed. It is the scope of this letter 
to present its ingredients and results of the simulations relevant for the pion 
form factor extraction via the radiative return method. 

Some partial results concerning the electron-positron pair production con- 
tributions to the pion form factor measurements and tests of the code were 
presented in [10,11], while in this letter the amplitude describing the reaction 
e + e~ — > e + e~7r + 7r~ is discussed, generation procedure is described in detail 
and results based on the complete tree level amplitude are presented. 



2 The scattering amplitude and the generation procedure 

As stated already in the introduction, the reaction e + e~ — > e + e~7r + 7r~ plays 
a role in the pion form factor measurement via the radiative return method 
only if the photons in the reaction e + e _ — > 7t + tt~ + photons are not tagged. 
This version of the radiative return method was used already by KLOE [2] 
and as more accurate analysis, based on a significantly bigger data sample, 
is expected, a detailed study of all possible contributions is obligatory. The 
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complete set of the lowest order diagrams, describing the reaction e + e 
e + e~7i + TT~ , is shown schematically in Fig.l. 








Fig. 1. Diagrams contributing to the process 

e + {pi)e~{p2) — > vr + (7ri)7r~(7T2)e + ((?i)e~(g2): initial state pair emission (a), fi- 
nal state pair emission (b,c), pion pair emission from t-channel Bhabha process 
(d) and 7*7* pion pair production (e). Only one representative diagram for a given 
set of diagrams is shown. 

Helicity amplitude method, with the conventions described in [12,13], was 
used for the scattering amplitude evaluation. It allows for a fast numerical 
evaluation and, in addition, all interferences are easily included. Moreover, it 
partly avoids numerical cancellations present, when one uses the trace method 
to get the square of the amplitude. To model photon-pion interactions, we 
use scalar QED (sQED) combined with the vector dominance model (VDM). 
Within these assumptions, the amplitude has the form 



M = M a + M b + M c + M d + M e 



(2) 



where the amplitudes Mj 



, e correspond to the contributions from 



the diagrams (a)-(e) from Fig.l. They read 
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«e 4 , J(VQ-2p 1 -T)Y ^{VQ-2 qi -T)\ . . , . 



tiQ 2 ' 7 V - 2 Q ■ P2 Q 2 + 2Q-q 2 



M e = ^F(t)F( tl ) ^ {pi)Yv(qiMq2hM (7) 

| 2v{p 1 )^ 1 v(q 1 )u(q 2 )^2u{p2) 2v(p 1 )f 2 v(q 1 )u(q 2 )f 1 u(p 2 ) \ 
*i - 27Ti(pi - gi) t + 27r!(g 2 - p 2 ) / 

where 

= F(g 2 )(7n - 7T 2 f , fci = gi + g 2 , Q = m + vr 2 , (8) 
s = (pi+ P2) 2 , * = (g 2 - P2) 2 , h = (pi - gi) 2 (9) 

and F(Q 2 ) is the pion form factor, which was adopted from [14]. For the 
numerical evaluation, following again the method described in papers [12,13], 
the amplitude is rewritten into a form where only 2x2 matrices appear. 

We use the multi-channel variance reduction method to improve efficiency of 
the generator and the generation is split into four channels, where two of them 
absorb peaks present in t-channel diagrams and other two take care of the 
s-channel peaks. 

For the t-channel peaks absorption, we use the following phase space repre- 
sentation 

dLips 4 (pi + p 2 ;q 1 ,q 2 ,7T 1 , tt 2 ) = 

dQ' 2 dQ 2 
dLips 2 (pi + p 2 ; Q', q 2 )— — dLips 2 (<2'; Q, gi)-— dLips 2 (Q; m, tt 2 ) (10) 

in one of the channels and an analogous one, with q 1 <-> q 2 , in the other 
channel. As both channels are completely symmetric under qi <-> q 2 , we will 
describe here only changes of variables, which smoothen the distribution, only 
in one of them. For the two invariant masses (Q 2 and Q' 2 ) the following change 
of variables was performed 
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Q 2 = (Vs - 2m e - v 73 ^) , z = --( v G-2m e -2m n ) 2 (l-r Q 2) , (11) 



CT = ^ + ml, (12) 
2/ = ~3(g 2 + 2 v ^m e ) 3 + (3(g 2 + 2v^m e ) 3 ~~ 3(s - 2^^) rQ ' 2 ' 



The angles of g 2 vector are defined in the initial e + e~ center of mass (cms) 
frame with z-axis along pi and the polar angle is used to absorb the peak 
coming from the propagator of the photon exchanged in the t-channel 



3m 2 . - s + Q' 2 - 2t ^ 1 
cos ^ 2 = / j=, — : — - , t = 



1 | sy/l-*£\V*(s,<p ,ml) 
t- ml(Q' 2 - m 2 e y 



y=~— + ZTTTT^ ZT^ r d n , K = 27ir^ M , (13) 



where 



t. = \ Uml-s + Q' 2 -^l-^\^(s,Q' 2 ,ml)) . (14) 



The angles of the Q vector are defined in Q' rest frame and the appropriate 
change of variables reads 



2E[Q -Q 2 - w jQ l3 

2|p1||0| 



^ Q = , <P Q = 2vrr, q (15) 



-1 2 

X ~ 2^1101(^-2^0-2^1101) + (Q 2 -2^g ) 2 -4b' 1 | 2 |0| 2reQ ' 



wnere^- 2 ^ , ^ - 2 ^ , - , M- 2 ^ 

and A(a, 6, c) = a 2 + 6 2 + c 2 — 2(ab + ac + be). As we chose here the z-axis 
along the p[ (the pi in the Q' rest frame) the vectors have to be rotated after 
generation to restore the general choice of the z-axis along pi in the e + e~ cms 
frame. 

Finally the angles of the positively charged pion are generated in the Q rest 
frame with flat distributions 
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cos 9 ni = -1 + 2r 07Ti , <p Wl = . 



(16) 



The described change of variables transforms the phase space into a unit 
hypercube (0 < rj < 1 , i — Q 2 , ■ ■ ■ , (f> ni ) and collecting all the jacobians 
it reads 



dLips 4 (pi + p 2 ; qi, q 2 , n u n 2 ) = P{q u q 2 )dr Q 2dr QI 2dre q2 dr 4)q2 dre Q dr 4)Q dre^dr 4> ^ 

(17) 

with 



Ami 



t 2 (Q' 2 - m 2 e ) 2 (Q 2 - 2Q ■ Pl f VQ^(V~s - 2m e - 2m n ) 2 



(Q 2 - 2E[Q y - 4|K| 2 |Q| 2 v^-H-v^ 
1 1 \ 



(g 2 + 2v^m e )3 {s-2^m e y 



(18) 



For the s-channel generation it is convenient to write the phase space in the 
following form 



dLips 4 (pi + p 2 ; qi, q 2 , 7T 2 ) = 

dk 2 dQ 2 
dLips 2 (pi + p 2 ; Q, fci)-^dLips 2 (A;i; q u g 2 )— dLips 2 (Q; m, tt 2 ). (19) 



The two generation channels used here differ only in the generation of the 
electron-positron pair invariant mass k\ and the change of variables will be 
described simultaneously. The invariant mass Q 2 is generated with a flat dis- 
tribution 



Q 2 = Ami + ({Vs - Zm e ) 2 - 4m^)r Q2 . 



(20) 



Reflecting two leading k\ behaviours of the cross section, two distinct changes 
of variables are done in the generation of k\: 



l/3x 



k 1 =sexp(y I 
y I = ln 3 (Am 2 Js) + 



(21) 



\n 3 ttl - y/cpjtf) -ln 3 (Am 2 Js) 



r kli 
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ki = s(l-eM-yn)) , (22) 

/v / Q 7 (2 v / i-v / Q 7 )\ 
VII = - ln(l - Amjs) - In I (g _ I r fc 2 j7/ . 



The &i polar angle is used to absorb peaks coming from the electron propa- 
gator, while its azimuthal angle is generated with a flat distribution: 



(p kl =27rr 0fci , cos# fcl = — 2|fci||pi| \2J 

, ( kj-2k w p w + 2\k 1 \\p 1 \ \ / ^-2/c 10 p 10 -2|fc 1 ||p 1 | \ 2 
^ m - 2k wPl0 - ^hM) + m [kj - 2k loPw + 21^1) n ^ [A6) 

where k w = *±g^, p 10 = & \p,\ = and \k\\ = are 

defined in the p\ + p 2 rest frame. 

The qi and the 7fi angles are generated with flat distributions 

4> qi = 2-kt^ , cos 9 qi = -1 + 2r 6qi , Wl = 27^ , cos ni = -1 + 27^ . (24) 

After the described changes of variables are performed, the phase space reads 
(i = I or //) 

dLips 4 (pi +p 2 ; 91,92,71-1,^2) = P s ^ i dr k 2dr Q 2dre ki dr^ ki dre qi dr 4>n dre^dr^ i ,(25) 
with 

p - = i^\/ 1 -|rv /l - ! i l ^.« 2 . *» (<VS - W - to") 



2fcioPio - ^i V 2 l^ilbi| "V V 2 |/ci||pi 

' -fc| + 2fciopio + 2|fci||pi 
-k\ + 2/c 10 p 10 - 2|fci||pi 



•P.-lnf-f + 2 ^o + 2 fi HP}h , (26) 



where 



« = ((i - MI) 2 ) - W / s) , ^ - m ( ^ffigj w 



for the change of variables from Eq.(21) or Eq.(22) respectively. Again < 



Ti < 1 for i = kf, 
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The function, which approximates the peaking behaviour of the matrix ele- 
ment reads 



— , with P. 



P s j + bP aJI 



31n 2 (fcf/ s ) 

ri r 



(28) 



The introduced a priori weights a and b, which guarantee the right relative 
contributions from different generation channels, were set to a = 1.1 and 
b = 1000, a choice optimal for DAPHNE energy. 



3 The pair production and the radiative return 



As the only experiment, which uses the radiative return method without pho- 
ton tagging is KLOE, we present here the results for DAPHNE energy only. 
The cross section of the reaction e + e~ — > 7r + 7r~ + 7(7) was obtained with the 
PHOKHARA 5.0 [8] event generator. For the event selection used by KLOE 
[2] the relative contribution of the pair production to the cross section of the 
reaction e + e~ — > 7r + 7r~ + 7(7) is shown in Fig. 2 (left). It amounts up to 
1.5% in the vicinity of the production threshold, but it is below 1% in the Q 2 
region, where the measurement was performed [2] (Q 2 > 0.33GeV 2 ). 
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Fig. 2. The ratio of the differential cross sections of the reactions e + e~ — > e + e~Tr + ir~ 
and e + e~ — > 7r + 7r _ 7(7) for KLOE [2] event selection (left) and separate contribu- 
tions to the cross section of the reaction e + e~ — * e + e~n + ir~: s - diagrams from 
Fig. la-lc included; s + t - diagrams from Fig. la-Id included; s + t + 27* - all 
diagrams from Fig. 1 included (right). 

For small Q 2 values the main contribution is given by the t— channel dia- 
grams, while around the p resonance the s-channel dominates (see Fig.2(right) 
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Fig. 3. The ratio of the t-channel and the s-channel contributions to the differen- 
tial cross section of the reaction e + e~ — ► 7r + 7r _ e + e _ (left); the relative contribu- 
tion of the photon fusion process to the differential cross section of the reaction 
e + e~ — * ir + ir~e + e~ (right) 



and Fig. 3). In the latter it was shown in [10], the ISR contributions 

dominate and the FSR is irrelevant. Moreover, the contributions to the cross 
sections from diagrams from Fig.lc are always small and the two-photon con- 
tributions (Fig.le) are completely negligible for the KLOE event selection as 
already pointed out in [15] and shown in Fig. 3 (right). This may not be true 
for different event selections, as the photon fusion diagrams can give sizable 
contributions especially in low Q 2 region. An example is presented in Fig.4, 
where only cuts on pion angles are applied. The pair production cross sec- 
tion can be in this case almost as big as the e + e~ — > 7r + 7r~7(7) cross section, 
where the s + t channel diagrams (Fig.l(a-d)) give only up to 3% of the 
e + e~ — > 7r + 7r~7(7) cross section [11] and the main contribution is given by 
the photon fusion diagrams. 



In the case of s-channel ISR and the pion pair emission from t-channel Bhabha 
diagrams a factorization analogous to the one of Eq.(l) occurs, with the ra- 
diator function given by QED. As a result, for the case of the KLOE event 
selection without photon tagging, it is possible to use the radiative return 
method adding contributions from the photon(s) and the pair production. 
The pion form factor, to be extracted from the data, is the same in both cases 
and the radiator function is a sum of both contributions. This procedure will 
be necessary, when the accuracy of the measurement will be below 1% at low 
Q 2 values. Alternatively, one can treat the pair production as a backgroud and 
study carefully the accuracy of its estimation. 
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Fig. 4. Various contributions to differential cross sections of the reaction 
e + e~ — ► e + e _ 7r + 7r _ : s - diagrams from Fig. la-lc included; s + t — diagrams 
from Fig. la-Id included; s + t + 27* - all diagrams from Fig. 1 included (left); 
The ratio of the differential cross sections of the reactions e + e~ — > e + e~ir + ir~ and 
e + e~ — > 7r + 7r^7(7) (right). In both cases only cuts on pions polar angles are im- 
posed. 

4 Conclusions 



Basing on the developed Monte Carlo program EKHARA presented in this 
letter, it was shown that the reaction e + e _ — > e + e~7r + 7r - may give non neg- 
ligible contributions to the pion form factor measurement via the radiative 
return method without photon tagging. For low invariant masses of the two- 
pion system it is up to 1% for the event selection used in KLOE analysis [2], 
but it can be substantially larger if some of the cuts, used in the analysis, are 
relaxed. 
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